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Abstract. Black hole initial data is usually produced using Bowen-York type puncture initial data or by applying an excision 
boundary condition. The benefits of the Bowen-York initial data are the ability to specify the spin and momentum of the system as 
parameters of the initial data. In an attempt to extend these benefits to other formulations of the Einstein constraints, the puncture 
method is reformulated using distributions as source terms. It is shown how the Bowen-York puncture black hole initial data and 
the trumpet variation is generated by distributional sources. A heuristic argument is presented to argue that these sources are the 
general sources of spin and momentum. In order to clarify the meaning of other distributional sources, an exact family of initial 
data with generalized sources to the Hamiltonian constraint are studied; spinning trumpet black hole initial data and black hole 
initial data with higher order momentum sources are also studied. 



1. Introduction 



The necessity that gravitational collapse leads to the generic formation of singularities distinguishes general relativity 
from other classical theories where point sources only enter as idealizations. Unfortunately, the use of distributions 
to describe sources in Einstein's equati ons is limited to a subset of distributions which is just regular enough for the 
equations to be well defined as shown in lGeroch and TraschenI (Il987b . These distributions must not be too singular and 
have at most codimension 1 . Nevertheless, many studies have described the structure o f various singular sp acetimes using 
linear distribution theory or generalized function algebras; these studies are reviewed in lSteinbauer and Vickers.(.2006) . In 
particular. IParkeii (1 1 979h characterized the stress energy of the Schwarzschild spacetime in Kerr-Schild coordinates using 
Schwartz distributions while the Kerr solution has been characterized usin g Colombeau algebras(a differential algebra 
which seeks to provide a rigorous multiplication of generalized functions) in |Parkeii(ll979h . 

In spite of the difficulties involved with the use of distributions in general relativity, distributions have been implicitly 
used in the construction of initial data for quite some time; this will be shown in detail in this paper. In particular, the 
puncture family of initial data (including the trumpet type initial data) can be interpreted as arising from differential 
equations with distributional sources. The usual derivation of puncture initial data considers that the data exists only on a 
punctured plane M."^ — {^i} with {xj} the locations of an arbit rary numb er of punctures; these locations are interpretable 
as an asymptotically flat extended sheet of the spacetime (see iBrandt and Briigmann ( 1997) ) and therefore as coordinate 
singularities. Boundary conditions are not imposed at these inner boundaries, which leads to an essential non-uniqueness 
to the constraint equations: considering the linear part of the differential operator on the entire plane R'^ with distributional 
source yields functions laying in the kernel of the differential operator in the punctured plane. This is ordinarily remedied 
by implicitly applying a boundary condition; the conformal factor is decomposed in t o a siri gular part and a regular part, and 
the Hamiltonian constraint is extended to the entire plane JBrandt and Briigmannl (1199 7)). thereby imposing an implicit 
regularity boundary condition on the remainder while the singular part is chosen arbitrarily to be one of these distributional 
solutions to the linear part of the differential operator. 

To the extent that this "puncture" procedure is rigorous, problems of multiplication of distributions can be ignored 
and the distributions can be used to parameterise the functions lying in the kernel of the differential operator The resulting 
functions are regular on the punctured plane and solving for the remainder on the entire plane implicitly chooses a 
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boundary condition by regularity. In this paper, I follow this procedure to produce and study initial data with distributional 
"source" terms. 

The layout of the paper is as follows. In section|2]the Green's function for the momentum equation is derived in the 
conformally flat, maximally sliced transverse-traceless decomposition of the initial value problem for general relativity. It 
is also shown how the puncture solution to the Hamiltonian constraint arises from a distributional source. In section l23] the 
sources for the trumpet black hole initial data are derived. In section[3]exact initial data having dipole and quadrupole type 
mass moments are considered while in section|4]spinning trumpet initial data and intial data with higher order momentum 
moments are considered. Finally, in section |5] the results of the paper are summarised. Througout this paper, I use units 
where c = G = 1 and indices run over spatial indices unless otherwise specified. 

2. Solutions to the constraint equations parameterised by distributions 

2.1. Green's function for the momentum equation 

The Hamiltonian and momentum constraint of Einstein's equations in the ADM formulation are 

R + K^ - K,,K'' = 167rp (1) 

Vj{K'^ -g'^K) =8TTf (2) 

with R the Ricci scalar of the initial spatial hypersurface with metric gij and Kij it's extrinsic curvature with trace K. 
When one performs the conformal transverse-traceless decomposition constraint equations and chooses a conformally 
flat background gij = tp'^Sij (ip is the conformal factor) and maximal slicing K = 0, the Hamiltonian and momentum 
constraint are respectively 

V2V + 8"V"^^y^*^+27r/5 = 0, (3) 

V^V, + S^^d.djV^ - Srrj, = 0, (4) 

where p and ji are the rescaled energy density and momentum vector respectively given in terms of the normal to the 
initial hypersurface n" and the projection tensor 'jab and the stress energy tensor as 

p = 7/,5p = yj^Tabu'^n'' (5) 

h = ^'j« = ^'^Tafcn°7^ (6) 

(in these equations, sums are carried out over the spacetime indices) and Aij — djVi+diVj —2 / iSijdkV^ is the conformal 
extrinsic curvature with Kij = ■ip~'^Aij. Initial data satisfying these equations are commonly referred to as Bowen-York 
initial data. Since the conformal background is flat, indices are raised and lowered with the flat metric and in particular 
Vi = V^ . Because of the use of the flat background, coordinate delta functions will be used in the following analysis. 

The basic form of the stress energy tensor of point-like particles as an idealization of extended objects has been 
derived in a different context in ISteinhoff and Puetzfeldl (120101) . where neglecting higher order multipole moments, the 
stress energy tensor of a pointlike source was shown to be 

^/3^yab ^ u(ap6)j _ 5«=(ayb)v,5 (7) 

where the factor ^J—g is appearing because of the use of coordinate delta functions. By finding the fundamental solution 
to the momentum constraint, the source terms for Bowen-York type initial data will be derived and shown to be exactly 
analogous to these terms. 

The Malgrange-Ehrenpreis theorem guarantees the existence of a Green's function for the momentum constraint; to 
find it, the source is set to ji{x) — PiS{x) with its Fourier transform being ji = Pi. The fundamental solution is then 
found by Fourier transforming the momentum constraint, algebraically rearranging terms and integrating to get 
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This is of course the linear momentum solution first derived in lBowen and YorkI d 19801) showing that the source of linear 



momentum arises as a point source Dirac delta distribution. This result can immediately be used to write down the general 
solution to the momentum constraint (with homogeneous boundary conditions at infinity) 

This general solution can be used to show that the source j* = —l/2e^^'^SkdjS{x) with e*^*^ the flat space Levi-Civita 
tensor gives the vector field 

V. = ^^^, (10) 

which is precisely the angular momentum solution of Bowen and York. 

The combination of the delta function source for the linear momentum term and the gradient term for the angular 
momentum yields the total source for the momentum equation for puncture initial data 

J^=(p^'^e,,kS''^^^S{x). (11) 



Bowen and YorkI ([1980) noted that the parameters P' and 5* are exactly the ADM linear and angular momentum 
respectively. It is exactly this property that makes the Bowen- York initial data so useful. Heuristically, this can be expected 
to always be the case on a maximal slice, since the ADM linear and angular momenta (associated with the translational 
and rotational killing vectors S^"-) as limits of integrals on spheres Sr can be converted to volume integrals 

SttP, = lim / K^'b^'^dAa = lim / \/aiK''bS.'')dV (12) 

if one assumes that the extrinsic curvature is valid on the entire open ball B,-. The momentum equation (|4|i can then be 
used to simplify the volume integrals 



Pe = lim / jbCdV = / jbCdV = / jbCdV (13) 

''-^°o Jb, Jb,. Jb^ 

in the second equality the limit vanishes because of the compact support of ji, dV is the volume form in the initial 
hypersurface and dV being the volume form in the conformal background. When the only source term is a momenta, 
the translational killing vector measures precisely P* and when the only source term is a spin the rotational killing vector 
measures precisely S\ This argument suffers from the weakness of distribution theory in converting from the surface 
integral to the volume integral and from going from the invariant volume element to the coo rdinate one. It is com pletely 



valid however within the nonlinear generalized function theory of Colombeau algebras (see [Grosser et alJ (1200 Ih for an 
introduction to this theory). 

2.2. The Hamiltonian constraint 

The usual way for solving the Hamiltonian constraint (O in the puncture formalism is by decomposing -0 = '(/'singular + u 
and having each component satisfy the following equations respectively 

V^Vsingular + 27rp = 0, (14) 

V^u + 8-'^ip-''A,jA'^ = 0, (15) 

where (fT4b is implicitly solved using a pointlike source p — MS(x — xj) so that '(/'singular = 1 + M/2\x — xi\. Then 
one needs to solve only the single e quation (usuall y numerically) for puncture initial data, ( fTSl ). For vanishing extrinsic 



curvature, it is already known from Parkeri (1197 9) that the source term for Schwarzschild solution is a delta function, 
although there it is associated with a physical singularity while here it is associated with a coordinate singularity. 

One should note that the original problem is ill posed (the solution is not unique) on the punctured plane M.^ — {xj} 
without suitable boundary conditions at the origin. Namely, any solution to ( fT4b having distributional source p at the 
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punctures lies in the kernel of the Laplacian; for instance the mass term M above is not unique-any mass will solve the 
Hamiltonian constraint ([3]). This ill-posedness could be removed by considering the original solution on the entire plane 
and considering the confomal factor to have a Dirac delta source term, in this case the puncture formalism is just a form 
of regularization but a nonlinear theory of distributions is necessary. Alternately, the above decomposition and implied 
regularity on u can be considered in lieu of a boundary condition. This ill -posedness is oft en exploited to iterate over 
the bare mass until the black hole has a desired horizon mass as outlined in lBriigmann et al.l (120081) . Parameterising the 
singular part via a distribution is preferable to appealing to a topological condition such as a compactified alternate world 
sheet if only for its simplicity. 



2.3. Trumpet initial data 



It was numerically verified in iHannam et alJ (120071) that the evolution under dynamic gauge conditions of puncture 
spacetimes causes the power of the singularity of the conformal factor to change in time. For this reason, the exact slicing 
of the Schwar zschild geometry that sat i sfies a gauge condition close to the normal gauge condition use d in evolution 
was found by Baumgarte and Naculich (2007). This solution has been used to produce initial data in IHannam et al. 



(2009); I mmerman and Baumgarte (2009) on the punctured plane which has been termed "trumpet" initial data. It should 
be emphasized that with zero momentum the trumpet black hole is identically a Schwarzschild black hole. Since the 
formulation used is again Bowen and York's conformally flat, transverse traceless decomposition, the singular structure 
which is implictly parameterising the initial data can be studied using the Green's function from the previous section. 

For a trumpet black hole, there is a background ext rinsic curvature in addition to the B owen- York spin and momentum 
terms. This background extrinsic curvature is given bv llmmerman and Baumgartd (l2009h 



Aij — 



4^3 



(^». - ^) 



and the corresponding exact conformal factor is 



^ 



R{t) 



with the areal radius R{r) given implicitly in terms of the isotropic radius r 



2R + M + V4i?2 + AMR + 3AP 



(4 + 3^/2)(2i?-3Af) 



8R + MI + 3VSR^ + SMR + 6M^ 



l/\/2 



(16) 



(17) 



(18) 



The presence of terms in the extrinsic curvature which scale inversely as r^ is an indication that the momentum 
sources must contain derivatives of the delta function. A generic source of the form 

j, = Bld,S{x) (19) 

together with the general solution to the momentum constraint (O can be used to determine the corresponding first order 
solution to the momentum constraint. Coefficients of a;"2/™z^ are collected to determine the unknown coefficients Bf 
such that (fT6l l is reproduced by this source term. In this way, the trumpet momentum source is determined to be 

VsaP 



Jl 



-d,5{x - xi) 



Interestingly, the trumpet initial data has no source for the Hamiltonian constraint, 
constraint against some test function tp, one has. 



- / 2TTpLpdV = 



?AVV- 






<P 



dV 



(20) 
Integrating the Hamiltonian 

(21) 



where (yS is a compactly supported test function with support at the location of the singularity of -tp and dV is the flat 
volume element of the background metric. This integral is converted to a limit 



2TTpipdV = lim 



-B,. 



V'V^(^- 






■ys 



dV 



(22) 
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and integrated by parts 



- f 2TTpipdV = lim / i-dii^d'ip + '^/ y^l dV+ lim / il^'^dS (23) 

with dS the area element over the surface of the open ball Sr- Since the radial derivative of the test function is uniformly 
bounded, the surface integral is limited as 

dip 



lim d> i^^dS 



^o.Tsr- dr 



< sup 



•^-^O '5, 



lim 6 ipdS = (24) 



dr 
with the equality arising from computation with the conformal factor One more integration by parts yields 



2, , AjA'^\ „.^ ,. / di; 



- f27rpipdV = lim/ V^ + '^- U^dF- lim (/) (^--dS" = 0. (25) 

where the volume integral vanishes due to the momentum constraint (there are no sources outside of the open ball) and 
the surface integral vanishes by an analogous uniform bounding as above. 

In the case of the trumpet data, the source of the mass of the spacetime arises entirely out of an isotropic "momentum" 
source. For this reason, separating the Hamiltonian constraint to factor out the singular part of the conformal factor is 
necessary for trumpet data only to formulate a numerical method that does not have a singular solution whereas for the 
original puncture data it is necessary to deal with a fundamental ill posedness and select a mass term. In this case, where 
one factors out the singularity only for numerical purposes, it suffices to factor out only the singularity y%lM/2r in order 
to ensure that the remainder u is bounded (although a more sophisticated factor is needed for spinning trumpets), and the 
Hamiltonian constraint is well defined over all M.^. This fact is used in section |4] to produce spinning trumpet initial data 
as well as more general initial data. 

3. Initial data with higher order mass moments 

Far from a source, the source itself may be idealized as a point particle with various mass moments (usually described in 
terms of multipoles). In this case the only equation to be solved is the linear part of the Hamiltonian constraint (O. The 
matter distribution that enters the field equations in this case is a series of derivatives of the delta function. In Newtonian 
physics, where mass is always positive, a dipole moment (the first order moment) is always zero except when the first order 
moment is computed about an origin which is not the barycentre; it is unclear to what extent this is meaningful in general 
relativity. A quadrupole moment (the trace-free second order moment) quantifies the asymmetry of the distribution. 
Considering the first two matter moments, the conformal source term including up to quadrupole terms is 

p = (1 - D'd^ + R'^didj) 6{x) (26) 

with D' quantifying the dipole moment and R^^ quantifying the second order moment (not the trace-free quadrupole 
moment). For a generic matter distribution with unit mass, the terms in this idealization can be computed as D' == J px^dv 
and R^^ = J px^x^ dv. It is possible to derive many forms for the second order moment, however for a ring with negligible 
cross section 

R'' = liS'Ji + S161) (27) 

with scalar quadrupole moment q being a combination of the mass and radius of the ring. If the dipole moment is aligned 
to the z axis so that Z?' = (0, 0, d), the solution to the Hamiltonian constraint (O with source (|26] | is given by 

1 dz 



^ = l + 7^ + :^ + ^(^'+y'-2^') (28) 



so long as ji — and so Aij — Kij = 0. 

These intitial data generically contain coordinate singularities with spherical topology: a closed surface where the 
conformal factor goes to zero. In the physical space, these surfaces would correspond to extended points; all points on the 
coordinate surface are separated by zero distance. As seen in figure[T]and figure|2]these coordinate singularity surfaces are 
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generically contained within the apparent horizons of the spacetimes (see iThornburd d 1 9961) for an overview of apparent 
horizons and horizon finding). The interior regions of these coordinate singularities are topologically S^ . Perhaps more 
interestingly, these surfaces of equidistant points ensure that the origin is a finite distance away from any other point in the 
plane; much like the trumpet initial data, it is no longer possible to interpret the origin to be the conformally compactified 
infinity of an alternate world sheet. 
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Figure 1. Apparent horizons (outer curves) and coordinate singularities (inner curves) plotted for spacetimes with various dipole 
moments d added. The inner curves represent coordinate singularities in the metric where the spatial metric vanishes. The 
centroid of the apparent horizon shifts upwards with increasing dipole moment while the centroid of the coordinate singularity 
shifts downward. 



The qualitative effect of adding the multipole moments is to change the apparent horizon in a corresponding fashion. 
The dipole moment, for small values of d shifts the horizon in the respective direction while maintaining an overall 
spherical shape, as seen in figure[T] For larger values of the dipole moment, the apparent horizon becomes dimpled in the 
direction —I?*. When the value of the quadrupole moment is small, the overall shape of the apparent horizon is elliptical, 
as seen in the left plot of figure |2] For large magnitude quadrupole moments, the apparent horizon becomes indented on 
the equator for negative moments, and at the poles for positive moments. 

These deformations induced on the horizons by the added structure of the distributional source terms cause the 
horizon mass of the black hole to decrease with increasing magnitude of the moment. The mass of the black hole is 
identified with its irreducible mass since there is no angular momentum. 



A/i„. = 



A 



(29) 



The masses of the black holes considered in this section are plotted in figure [3] The difference between irreducible mass 
M\„ and ADM mass is an estimate for the total amount of radiation in the spacetime. From the figure, one sees that for the 
dipole black holes, where the horizon could be found, the maximum amount of radiation is less than 0.1%. The horizon 
mass of the black holes with quadrupole moments is symmetric with respect to the sign of the quadrupole moment, which 
is unexpected given the difference in the structure of the apparent horizon for positive and negative quadrupole moments. 
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Figure 2. Apparent horizons (outer curves) and coordinate singularities (inner curves) plotted for spacetimes with various 
quadrupole moments q added. For small deformations, as seen in the left plot, the horizons are nearly elliptical. For larger 
deformations, as seen in the right plot, indentations occur on the horizon. The inner curves represent coordinate singularities in 
the metric where the spatial metric vanishes. 



For the horizons found, the maximum radiation content is about 0.2% and 0.5% for positive and negative quadrupole 
moments respectively. 



4. Fully non-linear initial data with distributional sources 



Considering completely general sources for the extrinsic curvature necessitates solving the fully non-linear Hamiltonian 
equation. Considering derivatives of the Dirac delta to higher than second order will lead to stronger singularities 
appearing in the conformal factor-the m/2r term is insufficient to ensure regularity. This also occurs if one tries to 
solve the Hamiltonian equation for the bare trumpet conformal factor numerically. Furthermore, for spinning trumpet 
initial data, the exact implicit conformal factor is not strong enough to make the solution fully regular since the spin term 
(an antisymmetric derivative) enters at the same order as the mass term (a gradient). For numerical regularity for general 
sources, the following decomposition of the conformal factor is proposed: 



■0 = 1 



-^'' + ^ 



(30) 



Here tpp satisfies the Poisson equation V^V'p-|-27r/5 = 0, the term a depends only on angle, and u is the bounded remainder 
part of the conformal factor. The extrinsic curvature term AijA^^ can always be written in terms of a denominator which 
is a power of radius r and a numerator which is polynomial in r and trigonometric functions of the angles. Considering 
only axisymmetric sources and using spherical coordinates, the extrinsic curvature is decomposed as 



AijA 



K{r,9) 



Er=i^M^)+/?W 



(31) 
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Figure 3. Irreducible mass of non-spinning, zero momentum black hole with added dipole and quadrupole moments. The 
difference between the irreducible mass and the ADM mass A/adm = 1 estimates the fractional radiation content of these 
spacetimes. Spacetimes with dipole moments have very little radiation content in comparison to the spacetimes with quadrupole 
moments. The extent of the curves shows the range for which apparent horizons could still be numerically found numerically. 



with n the problem dependent integer such that /3 ^ 0. By demanding that lini r^^^V^w = 0, one concludes that 



2 — n = and the Hamiltonian constraint yields the following two equations 



Ao(7 + p(p- l)o- 



P 



lim r^ib, 

r-i-O ^ 



(32) 



r^+^y^u - 



lim r^ibo 

r-i-O '^ 



8 (rP (1 + u + Vp) + cr)^ 



= 



(33) 



with Aq — -§gj + cot 9-§g the angular part of the Laplacian multiplied by the radius squared. These equations would be 
analagous for non axisymmetric sources except that A^ would be replaced by the full angular part of the Laplacian. To 
study solutions to these equations, Chebyshev pseudo-spectral methods are used. In particular, the collocation points are 
chosen as the roots of the Chebyshev polynomials in order to ensure regularity at the origin and symmetry axis without 
the use of boundary conditions. In spherical coordinate s, the ordinar i ly C^ puncture remainder is much more regular 
and exponential convergence can be achieved as noted in lAnsorg et al.l (120041) . The exponential convergence is shown in 
figure m where the truncation error of a spinning puncture dataset displays an exponential index of convergence. 



4.1. Spinning trumpet initial data 

In the case of non-spinning trumpet initial data (with or without momentum) /3 = 81/8M^ which leads p = 1/2 and 
a = y%3A//2 so that the singular part of the conformal factor is V's = 1 + \/3M/2r. For spinning trumpets, it is 
insufficient to subtract this factor or even the exact trumpet conformal factor as shown in llmmerman and Baumgarte 
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Figure 4. The truncation en'or for a spinning puncture solution in splieiical coordinates with J = 0.5 and M = 1. The 
exponential fall-off with linear increase in basis order is indicative of the function having an exponential index of convergence. 
The use of spherical coordinates is crucial to achieve this result. 



(12009 1). However, it is possible to produce this initial data with the above decomposition where the coefficient of the 
singular part a has angular dependence. In the case of spinning trumpet initial data the extrinsic curvature is generated 
from the general solution (|9|l using the source 



= s^ 



Z'^IJZ 



d^ 



-d, 5{x). 



(34) 



Due to the fact that the limiting case of both trumpet and puncture initial data with spins will have the same source term 
ji = Sztijzdj5{x) the limiting behaviour of the solutions are identical. However, the transition from the puncture type 
l/r singularity to the spinning trumpet type 1/^ singularity is singular for punctures while the transition is smooth for 
the trumpet initial data produced by the method here. The exact singular solution may even be computed. It is known for 
Bowen-York ty pe initial data that it is not possible to produce black holes with arbitrarily large ratios of spin to mass; this 
was detailed in Lovelace et al.l (l2008h . To quantify the behaviour of the spinning trumpets, the dimensionless spin of the 
black hole 

is monitored at the initial hypersurface using a Christodoulou-like mass 



M.4 ■ 

irr 



(36) 



Due to the axisymmetry of the system, the spin measured on the horizon is identical to the ADM spin which is identical 
to the parameter Sz used to construct the extrinsic curvature. During evolution, the spin measured on the horizon will 
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Figure 5. The dimensionless spin of the black hole measured on the horizon x and at infinity ej as a function of bare 
dimensionless spin. Although the black hole has an initially large value of the dimensionless spin, it accretes gravitational 
radiation during evolution which increases its mass and subsequently decreases its dimensionless spin. This process causes the 
black hole to relax, under evolution, to a dimensionless spin close to the ADM value. 



quickly relax to a lower equilibriu m value and an excellen t approximation of the relaxed value is the dimensionless spin 
as measured at infinity as argued bv lLovelace et alJ (l2008h 

5. 



£j 



A/A2 



(37) 



here Sz is a parameter of the problem but all of the masses must be measured. These quantities are displayed in figure |5] 
where the bare parameters have been varied by many orders of magnitude. Exactly like the spinning puncture initial 
data, it is not possible to produce black holes of arbitrary spin parameter. Although the total angular momentum of the 
solution can be set to an arbitrary value, the mass of the spacetime increases yielding an asymptotic maximum value 
of the dimensionless spin. One sees that the limiting value of the dimensionless spin is x ~ 0.984 which will relax 
upon evolution to approximately the asymptotic value of ej w 0.937. The data presented here has been run at multiple 
resolutions to ensure accuracy to one part in one million, the data presented are results run with the number of collocation 
points set to (n^, ng) ^ (80, 40). 



4.2. Higher order momentum moments 

In a close binary, it is likely that the black holes will locally have a Kerr like geometry at first order, but due to tidal effects 
will have induced moments on their horizons. It may be desirable to incorporate these effects in the source terms used to 
construct initial data (or in the shape of apparent horizon). For this reason, in this section black holes with higher order 
momentum moments are considered. In particular, the following source is considered for the momentum constraint 



]t 



12 z 



d' + Qze,,zd^dz + Qce^Jke''"'^m^^ Six) 



(38) 



Distributional sources for black hole initial data 



11 



0.995 



0.990- 



^ 0.985- 



0.980 



0.975 



Q^ quadrupole 
Q^ quadrupole 



0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 
Q, or Oc 



Figure 6. The change in mass for a spinning puncture as the momentum quadrupole moments are increased. The irreducible mass 
is seen to increase with the quadrupole moment, from this one can conclude that the fractional radiation content in the spacetime 
is decreasing. 



with the factor 1/4 chosen to have a high but not numerically problematic spin xo ~ 0.423. The two terms are respectively 
an axial derivative of the spin term (proportional to Qz) and the curl of the spin term (proportional to Qc), both give rise to 
axisymmetric systems angular dependence of a higher frequency than the spin term. The Hamiltonian constraint is chosen 
to have source p = 5{x) to make the solutions as regular as possible. Due to the increased irregularity of the solutions, 
{ur, ng) = (120, 50) collocation points were used to achieve convergence over the range of parameters studied here. 

The effects of the quadrupole moments on the mass and spin of the black hole are plotted respectively in figure|6]and 
figure |2] respectively. The Christodoulou-like mass is not plotted for these spacetimes, as the black holes are no longer 
even approximately Kerr-like. It is seen from figure|6]that the irreducible mass of the black hole increases with increasing 
quadrupole moment. Since the ir reducible mass can never increase over the value of the ADM mass (this is the Penrose 
inequality, proved bv lBravl (1200 ih ) during evolution, there must be less radiation in the spacetime. Because of the increase 
in the irreducible mass of the black hole, the dimensionless spin of the black hole must decrease. This effect is shown 
in figure |7] where a very steep decrease in the spin is seen with an increase in the momentum quadrupole. This effect 
could have consequences for spinning black holes in a binary system; if the black holes are tidally interacting, the induced 
quadrupole moment could prevent the black hole from having an extremal spin by increasing its irreducible mass and with 
it enforcing cosmic censorship. 



5. Discussion 



The above analysis revealed that the puncture (and trumpet) initial data had as its origin distributional source terms. Two 
distinct terms that gave rise to mass terms were found, a delta function source for the Hamiltonian constraint in the case 
of puncture initial data and an isotropic delta function gradient source for the momentum constraints in the trumpet initial 
data. The linear momentum of the puncture initial data, and subsequently the trumpet initial data, arises from a delta 
function source while the angular momentum arises from the contraction of an anti-symmetric (spin) tensor with the 
gradient of a delta function. 
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Figure 7. The change in dimensionless spin for a spinning puncture as the momentum quadrupole moments are increased. Due 
to the increase in mass as seen in figure|6]the dimensionless spin decreases. 



It is possible to extend the method (parameterisation by distributions of the vacuum solutions) to the conformal thin 
sandwich formulation or to use alternate conformal backgrounds. The basic approach being to parameterise solutions 
lying in the kernel of the linear part of the differential equations, then to decompose the conformal factor into a source, 
singular and regular part and finally to solve for the regularised remainder on the entire plane. This process was used to 
generate novel black hole initial data in section |4] Given the heuristic argument presented in section |2l the same source 
terms would need to be used to generate the extrinsic curvature for linear and angular momentum. It is possible to put this 
heuristic argument on rigorous footing using Colombeau algebras but one must then tackle the issue of uniqueness with 
respect to mollifier if using special algebras. 

It is also possible to extend these results to the production of binary black hole initial data. For binary data, especially 
binary trumpet data, it may be preferable to produce individual spinning, boosted trumpets using the methods above to 
high order and then subtract the combined remainders so that the singular parts are correct to next to leading order The 
resulting remainder should be much more regular in thi s way. The higher ord er moments used here may become important 
in a binary, where tidal effects will distort the horizon jPoisson and Vlasovi (12010.) ) and with it possibly the sources. 

The code used in the production of the nonlinear initial data studied here is freely available under the Gnu General 
Public License at https://github.com/SwampWalker/LeapingMonkey. Readers interested in extending these results are 
welcome and encouraged to build upon the code used here. 



Acknowledgments 

I would like to thank Juan Barranco, Barry Wardell, Ian Hinder and Abraham Harte for their useful discussions and 
insight. Partial support comes from the DFG Grant SFB/Transregio 7. 



References 

Ansorg M, Brugmann B and Tichy W 2004 Phys. Rev. D 70(6), 06401 1-+. 



REFERENCES 13 

Baumgarte T W and Naculich S G 2007 Phys. Rev. D 75(6), 067502-+. 

Bowen J M and York, Jr. J W 1980 Phys. Rev. D 21, 2047-2056. 

Brandt S and Briigmann B 1997 Physical Review Letters 78, 3606-3609. 

Bray H L 2001 Journal of Differential Geometry 59, 111-261 . 

Brugmann B, Gonzalez J A, Hannam M, Husa S, Sperhake U and Tichy W 2008 Phys. Rev D 77(2), 024027-+. 

Geroch R and Traschen J 1987 Phys. Rev D 36, 1017-1031. 

Grosser M, Kunzinger M, Oberguggenberger M and Steinbauer R 2001 Geometric Theory of Generalized Functions with 
Applications to General Relativity Kluwer Academic Publishers Dordrecht, NL. 

Hannam M, Husa S, Brugmann B, Gonzalez J A, Sperhake U and Murchadha N O 2007 Journal of Physics Conference 
Series 66{\),0\2QA1-+. 

Hannam M, Husa S and Murchadha N 6 2009 Phys. Rev D 80(12), 124007-+. 

Immerman J D and Baumgarte T W 2009 Phys. Rev D 80(6), 061501-+. 

Lovelace G, Owen R, Pfeiffer H P and Chu T 2008 Phys. Rev D 78(8), 0840 17-+. 

Parker P E 1979 Journal of Mathematical Physics 20, 1423-1426. 

Poisson E and Vlasov I 2010 Phys. Rev D 81(2), 024029-+. 

Steinbauer R and Vickers J A 2006 Classical and Quantum Gravity 23, 91-+. 

Steinhoff J and Puetzfeld D 2010 Phys. Rev D 81(4), 044019-+. 

Thornburg J 1996 Phys. Rev D 54, 4899-4918. 



